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Abstract 

The annihilation of weakly interacting massive particles can provide an important 
heat source for the first (Pop. Ill) stars, potentially leading to a new phase of stellar 
evolution known as a "Dark Star" . When dark matter (DM) capture via scattering off 
of baryons is included, the luminosity from DM annihilation may dominate over the 
luminosity due to fusion, depending on the DM density and scattering cross-section. 
The influx of DM due to capture may thus prolong the lifetime of the Dark Stars. 
Comparison of DM luminosity with the Eddington luminosity for the star may constrain 
the stellar mass of zero metallicity stars; in this case DM will uniquely determine the 
mass of the first stars. Alternatively, if sufficiently massive Pop. Ill stars are found, 
they might be used to bound dark matter properties. 



1 Introduction 

The first stars in the Universe mark the end of the cosmic dark ages, reionize the Universe, 
and provide the enriched gas required for later stellar generations. They may also be 
important as precursors to black holes that coalesce and power bright early quasars. The 
first stars are thought to form inside dark matter halos of mass 1O 5 M0-1O 6 M0 at redshifts 
z = 10 — 50 [21 |3j. These halos consist of 85% dark matter and 15% baryons in the form 
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of metal-free hydrogen and helium gas. Theoretical calculations indicate that the baryonic 
matter cools and collapses via molecular hydrogen cooling [HOE] into a single small 
protostar [7J at the center of the halo (for reviews see e.g. [HI l9j [TO]). 

In this paper we continue our previous work [I] examining the effect of dark matter (DM)the 
particles on the first stars. We focus on Weakly Interacting Massive Particles (WIMPs), 
which are the favorite dark matter candidate of many physicists because they automat- 
ically provide approximately the right amount of dark matter, i.e. 24% of the current 
energy density of the Universe. Many WIMP candidates are their own antiparticles, in 
which case they can annihilate with themselves in the early universe with weak interaction 
cross sections, leaving behind this relic density. Probably the best example of a WIMP 
is the neutralino, which in many models is the lightest supersymmetric particle. The 
neutralino, the supersymmetric partner of the W, Z, and Higgs bosons, has the required 
weak interaction cross section and mass ~ GeV - TeV to give the correct amount of dark 
matter and would play an important role in the first stars. For reviews of SUSY and other 
dark matter candidates see [251 123 E7J . 

This same annihilation process is the basis of the work we consider here. WIMP self- 
annihilation is relevant wherever the WIMP density is sufficiently high. Such regimes 
include the early Universe, in galactic halos today [11[ I12j . in the Sun [13] and Earth 
|14[ I15j. and in the first stars. As our canonical values, we will use the standard value for 
the annihilation cross section 

Mann = 3 x KT 26 cm 3 /s, (1) 

as this gives the right WIMP relic density today, as well as take m x = 100 GeV for our 
canonical value of the WIMP particle mass, but will also consider a broader range of WIMP 
masses (1 GeV-10 TeV) and cross-section. 

The interaction strengths and masses of the neutralino depend on a large number of 
model parameters. In the minimal supergravity model, experimental and observational 
bounds restrict m x to 50 GeV-2 TeV, while (av) lies within an order of magnitude of 
3 x 10 -26 cm 3 /sec (except at the low end of the mass range where it could be several orders 
of magnitude smaller) \17\ [TH] . Nonthermal particles can have annihilation cross-sections 
that are many orders of magnitude larger (e.g. |16j ) and would have even more drastic 
effects. Given the present state of the field there are many types of DM candidates which 
could apply [19]; the effects we find apply equally well to other WIMP candidates with 
comparable cross sections for self-annihilation and scattering off of nucleons. 

In this paper we consider the effects of WIMP annihilation on the first stars. In a previous 
paper (hereafter, Paper I) [1] two of us (together with P. Gondolo) considered the effects 

*If the lightest supersymmetric particle were the axino or gravitino, it would not have weak interaction 
annihilations and would not produce the phenomena described in this paper. 
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of dark matter annihilation on the formation of the first stars. We found that a crucial 
transition takes place when the gas number density of the collapsing protostar exceeds a 
critical value (10 13 cm~ 3 for a 100 GeV WIMP mass): at this point WIMP annihilation 
heating dominates over all cooling mechanisms and prevents the further collapse of the 
star. We suggested that the very first stellar objects might be "dark stars," a new phase of 
stellar evolution in which the DM - while only supplying 1% of the mass density - provides 
the power source for the star through DM annihilation. We are as yet uncertain of the 
lifetime of these theoretical objects. Once the DM contained inside the star runs out, the 
star could contract and heat up to the point where fusion becomes possible. However, as 
the star reaches nuclear density, there is an additional mechanism for repopulating the DM 
inside the star: capture of more DM from the ambient medium. This new source of DM 
can extend the lifetime of the Dark Star. Indeed the capture process can continue as long 
as there is enough ambient DM passing through the star, causing DM annihilation inside 
the star to continue, possibly even to today. 

Dark Stars, powered by DM annihilation, require three key ingredients, as shown in Paper 
I: 1) high dark matter density, 2) annihilation products remain trapped in the star, and 3) 
DM heating dominates over the other sources of heating. The first stars exist at the right 
place and at the right time to have the best chance of achieving the first criterion: they 
exist in the high density centers of dark matter haloes, and they form at high redshifts 
(density scales as (1 + z) 3 ). These densities are still not enough: Paper I showed that 
the DM density must be driven up still further in order for the DM heating to become 
important: as the baryons condense to form stars, they come to dominate the potential 
well and pull the DM in with them and drive up the density. The resultant DM densities 
were computed in Paper I by adiabatic contraction. Once the DM due to this effect runs 
out, the Dark Star phase might end. On the other hand, there is another effect that can 
repopulate the DM inside the stars and allow the Dark Star phase to continue: capture of 
more DM from the ambient medium. To get a significant amount of captured material into 
the star requires the enhanced densities due to adiabatic contraction. We will estimate the 
required ambient DM density for capture to be important, and find that it is somewhat 
lower than is required in the earlier protostellar phase. 

In this paper, we consider the effects of DM annihilation on early zero-metallicity (Popula- 
tion III) stars, once they do have fusion inside their cores. These stars live inside a reservoir 
of WIMPs; as the WIMPs move through the stars, some of the WIMPs are captured by 
the stars. The captured DM sinks to the center of the stars, where the DM can annihilate 
very efficiently. This has the effect of dramatically increasing the annihilation rate inside 
of a star, compared to DM annihilation without scattering. The annihilation provides a 
heat source for the stars, and we compare this DM annihilation luminosity to the fusion 
luminosity of the stars, as well as to the Eddington luminosity. Of course, the presence of 
DM (with high densities due to adiabatic contraction or capture) would already become 
important during (and seriously affect) the formation and earlier stages of these stars. 
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Here it is simply our intent to show that under some circumstances DM capture makes the 
importance of DM heating in Pop. Ill stars unavoidable. Thus, even if DM annihilation 
fails to stop a Pop. Ill star from forming (as was considered in pQ), a Dark Star powered 
by DM annihilation may exist at the high nuclear densities where fusion can take place. 
Again, we wish to note that the DM supplies less than a percent of the mass of the Dark 
Star and yet may be responsible for its luminosity. 

The two key uncertainties in this work are: (i) the scattering cross section must be at 
(or near) the experimental limit, and (ii) the ambient DM density must be high enough 
for capture to take place. Whereas the annihilation cross section is fixed to be close to 
Eq. (PQ), on theoretical grounds the scattering cross section can vary across many orders of 
magnitude. As discussed later, it is, however, constrained by experimental bounds. For 
scattering to matter in the first stars, the scattering cross section must be within two 
orders of magnitude of the current experimental limits. Such a cross section should be 
experimentally accessible in the next round of DM detection experiments. The second 
criterion is likely to be true for a while after the star is created, but it is not clear for how 
long. Once the halo containing the Dark Star merges with other objects, it is not clear how 
long the central DM in the halo remains undisturbed, and it is not clear how long the Dark 
Star remains at this central point. In principle the capture could continue indefinitely so 
that a dark star could still exist today, but this is very unclear. 

Previous work on DM annihilation powering stars has also been done in the context of high 
DM densities near the supermassive black holes in galactic centers, e.g. WIMP burners 
[20] and more generally [21] [22] . 

Just as we were preparing to submit our paper, a very similar work was submitted by [23J . 
We have carried the analysis further. We agree with [23] in the conclusion that DM 
annihilation may dominate over fusion, and we illustrate the DM densities required for 
this conclusion. In addition, we go one step further and discuss the possibility that the 
DM power source may exceed the Eddington luminosity and prevent the first stars from 
growing beyond a limited mass. This would effect the IR background, the re-ionization of 
the universe, the number of supernova, and potentially the nature of supernova of the first 
stars; this will be addressed in a separate publication [29] . 

We begin by discussing the equilibrium WIMP abundance in the first stars, by computing 
the number of WIMPs captured by the first stars, and equating this with the annihilation 
rate of WIMPs in the first stars. A discussion of adiabatic contraction, which may drive 
up the DM density near the baryons, follows. Then we compute the DM annihilation 
luminosity and compare with fusion luminosity. Finally, we compare the DM luminosity 
with the Eddington luminosity and find a maximum stellar mass as a function of DM 
density. 
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2 WIMP Abundance 



As WIMPs travel through a star, they can scatter off of the nuclei in the star with the 
scattering cross section a c . Although most of the WIMPs travel right through the star, 
some of them lose enough energy to be captured. We call the capture rate C(s~ 1 ). The 
WIMPs then sink to the center of the star, where they can annihilate with one another 
with the annihilation rate r^(s _1 ). This process was previously noticed as important for 
the Sun by [13] and in the Earth by [141 115] . DM could be detected through neutrinos 
produced from the DM annihilation products in the Sun or Earth with experiments such as 
SuperKamiokande [32] and AMANDA (which did not find a signal and placed bounds) and 
ICECUBE (which is starting to take data)[30j. Other indirect searches such as GLAST 
and PAMELA could detect the gamma-ray and positron annihilation products respectively 
from DM annihilating in the Milky Way halo (the gamma-rays and positrons from DM 
annihilating in the Sun or Earth would be trapped inside the objects and would not make 
it to our detectors). 

The number of WIMPs N in the star is then determined by a competition between capture 
and annihilation via the differential equation, 

N = C - 2F A = C - C A N 2 , (2) 

where T A is the annihilation rate (and the factor of two appears because two particles are 
annihilated in each event), and 

C A = 2T A /N 2 (3) 

is defined as an N— independent annihilation coefficient. Solving this equation, we find 
1^4 = 7jCtanh 2 (t/r), where 

r = (CCa)- 1 ' 2 (4) 

is the equilibration timescale. Equilibrium corresponds to a balance between the capture 
and annihilation rates, i.e. 

Ta = \c (5) 

As we will show, for the case of Pop III stars, equilibrium is quickly reached (t 3> r) and 
we may use Eq.([5]). 



2.1 Annihilation Rate 

The annihilation rate 



r A 



J d 3 rn x (r) 2 (av) ann (6) 
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where n x (r) is the density of captured DM at a point r inside of the star. In this section 
we will show that the WIMPs quickly thermalize with the core of the star, so that we can 
treat the DM density distribution as isothermal, and then we will compute the annihilation 
rate. Our work closely follows the approach previously given by [24J. 

First let us examine the WIMP thermalization timescale inside the star. The amount of 
energy AE lost by a WIMP in a scattering event with a nucleus (proton mass M p ) in the star 
ranges from < < ^ m m ^M P )/2} 12 • Assuming a flat distribution and taking M p <C m x , 
the average energy loss is 2M p /m x . Roughly, thermalization requires AE/E ~ 1; i.e. there 
must be m x /2M p scatters. Thus, for m x S> 1 GeV, the timescale for thermalization can 
be estimated as 

1 M x 

T th « — f- (7) 

a c v esc riH 2M H 

where v esc is the escape velocity of a DM particle from the surface of the star, and nn 
is the average density of the star. For a hundred GeV WIMP, a c = 10~ 39 cm~ 2 , and 
tih = 10 cm -3 , we find that the thermalization time scale to be very short, roughly three 
months. 

Thus we can use an isothermal distribution for the DM, 

n(r) x = n c e- m ^ /kT (8) 

where n c is the central number density of DM and T is the central temperature of the 
star, 

f r GM(r) , 

0(r) = / y —^dr. (9) 

Jo r 

is the gravitational potential at radius r with respect to the center, and M(r) is the mass 
interior to r. One can define effective volumes 

rR* 

Vj = 4tt / r 2 e- jm ^ T dr. (10) 



o 



Upon integration this gives 



Vj = [3m 2 pl T/(2jm x p c )f 2 , (11) 



where m p \ is the Planck mass, and p c is the core mass density of the star. The name 
"effective volume" is suggestive since we have iV = n V\ and also the total annihilation 
rate is given as T = {crv) !mn n 2 1 V2- 



One can then solve Eq.([2|) to find the N-independent annihilation coefficient defined in 
Eq.©, 



C A = Mann-TTl- (12) 
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M*(M ) 


T(K) 


/?(g/cm 3 ) 


V 2 /V?(cm 


" 3 ) 


Ca{s- v 


) 


Sun 






1.72 x 10" 


-28 


5.16 x 10- 


-54 


10 


9.55 x 10 7 


225.8 


1.77 x 1(T 


-29 


5.31 x 10- 


-55 


50 


1.13 x 10 s 


48.63 


1.38 x 10" 


-30 


4.14 x 10- 


-56 


100 


1.18 x 10 s 


31.88 


6.86 x 1(T 


-31 


2.06 x 10- 


-56 


250 


1.23 x 10 s 


19.72 


3.11 x 10- 


-31 


9.33 x 10- 


-57 



Table 1: Central temperature T(K) and central baryon density p(g/cm s ) for various 
masses of metal free stars half way through hydrogen burning |28j . The effective volume 
V2/V 2 is also shown. The DM annihilation rate is T a = ^C^iV 2 where N is the number 
of WIMPs in the star. Please note that the entry marked "Sun" refers to the present 
day Sun for comparison. We have used the fiducial values m x = 100 GeV and (cr?;) ann = 
3 x 10 _26 cm 3 /s as needed. 

In Table 1, we have computed Ca in Pop III stars of different masses. From S. Woosley [28] . 
we have obtained the properties of zero metallicity stars when they are halfway through 
hydrogen burning on the main sequence. In the table we have used our canonical values in 
Eq.([T]) for the annihilation cross section and mass (m x = 100 GeV); the results easily be 
scaled to other values since Ca oc m}; 5 {av) ann . 



2.2 Capture Rate 

WIMP interactions with nuclei are of two kinds: spin-independent, which scale as A 2 
(where A is the number of nucleons in the nucleus), and spin-dependent, which require the 
nucleon to have a spin. Currently the experimental bounds on elastic scattering are the 
weakest for the spin-dependent (SD) contribution (to be precise we are considering only 
scattering off of protons since the stars are comprised primarily of hydrogen). In figure 
2 [31] we illustrate bounds on the SD component from direct detection as well as from 
Super-Kamikande jSU [32j [39l [33J [35], EU [37] . The latter, which are the most constraining, 
assume that a significant fraction of the annihilation energy goes into neutrinos (as is likely 
for SUSY particles); if the neutrino component is small then the SD cross section could 
be several orders of magnitude higher. As our fiducial value, in this paper, we use the 
spin-dependent cross section 

a c = 10~ 39 cm 2 (13) 

which is consistent with all experimental bounds, but we will always show the dependence 
of any result on the value of a c . The bound on the spin-independent (SI) scattering is 
much tighter, a si < 10 _42 cm 2 for m x = 100 GeV |38|. 139]. The first stars are made only 
of hydrogen or helium, so the A 2 enhancement for the spin-independent contribution is 
not substantial. In this paper we consider only the spin-dependent contribution, though in 
principle for any specific candidate WIMP one should self-consistently include both. 
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The capture rate per unit volume at a distance r from the center of the star, for an observer 
at rest with respect to the WIMP distribution (as should be a good approximation here), 
is HUHO] 



dC , , /6\ 1/2 ... Mr) 2 



(r) = - n(r)n x {r)(a c v) 



ok \tt I v z 



1 



1 - exp(-B 2 ) 

S2 



(14) 



where n is the number density of nucleons (here, hydrogen), n x is the WIMP number 
density, v(r) is the escape velocity of WIMPs from the star at a given radius r, v 2 = 
is a "velocity dispersion" of WIMPs in the DM halo, and 

2 _ 3v(rf n 
B =2 — ~ (15) 



where \i = is the ratio of WIMP to nucleon mass and /i_ = (fx+l)/2. For an observer 
moving with respect to the WIMPs, the quantity in square brackets in Eq. (|14p becomes a 
more complicated function of B and the relative velocity as shown in Eq.(2.24) of |40j . 

The capture rate for the entire star is then 



" „ , dC 



C= \ 4irr 2 dr — (r), (16) 



o 



where i?* is the radius of the star. To obtain a conservative and fairly accurate estimate 
of the capture rat^l , we may take 

v(r) 2 = v(R+) 2 = 2 ^k^ vL (17) 

for all r, assume the term in square brackets in Eq. (jl4|) is very close to 1 (justified below), 
and take a uniform dark matter density. In this case the integral simplifies to give 

c=^y /2 (^/„w,(^) 2 ^, (18) 

where M* is the stellar mass, m p is the proton mass and fn is the fraction of the star 
in hydrogen. (Note that hydrogen has spin while helium generally does not. We could, 
in principle, consider the spin-independent contribution of scattering off of hydrogen and 
helium in the stars; we have not done so because we believe this contribution to be subdom- 
inant. In any case the current work is conservative in considering only the spin-dependent 
scattering.) 



2 We have subsequently obtained much more accurate values for v(r) by treating the Pop III star as 
a polytrope with index n = 3 (a good approximation), and confirmed that the conservative capture rate 
presented here differs (it is too low) by only a factor of a few. 
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To estimate v as per [53], we take the virial velocity of the DM halo, 

\W 



<C2> =I4±' (19) 



where 



W = -4irG J PhaioM ha i (r)rdr (20) 

and the typical DM halo containing a Pop III star has Mh a i = 10 5 — 1O 6 M0. 
We use a Navarro, Frenk & White (NFW) profile 06] for the DM, 

Phalo = — — — — ~2, (21) 

where r s is the scale radius. The normalization po, known as the central density, depends 
on the concentration parameter C v ; r and on the redshift when the Halo virializes Z v { r . 
These parameters range from C v ; r = (1 — 10) [32] and Z v - W = 10 — 50 [21 [3] . With r s in the 
range (15-100) pc, we find v = (1 — 15) km/s. As our fiducial value, we will take 

v = 10 km/s. (22) 



For B 1, the term in square brackets in Eq. (|14p is very close to 1, and this holds for 
all stellar and WIMP masses we are considering. (For example, for a 1M & star, we have 
^esc = 618 km/s, which is much larger than v = 10 km/s; then for m x = 100 GeV we 
find B ~ 100.) Thus we may ignore the term in square brackets. The bracketed term 
changes for a star moving through the WIMP halo (rather than being stationary as we 
have assumed), but the factor that replaces the term in square brackets is 0(1). (For 
example, for a 1 Mq star, with m x = 100 GeV, and with a velocity v, we find that the 
factor is 0.66.) We note that in today's stars this factor is much more important than in 
the first stars because v is much larger today (due to the fact that today's galactic haloes 
are much larger, e.g. 1O 12 M0). 

In Table 2, we evaluate the capture rate in Eq. (fT8|) . again using properties (including stellar 
radius) of Pop III stars from [28] . In obtaining these numbers, we have used p x = 10 9 
GeV/cm 3 , m x =100 GeV, and a c = 10 _39 cm 2 ; the result can easily be scaled to other 
values since C oc p x a c /m x . For general values we find that: 

C = 49xl0 3 V 1 ( Vcsc Y( * V 1 ( ^ ^ ( mx Y 1 ( ^ 1 

\M Q ) V618km/sy V10km/sy \ 10 9 GeV/cmy V 100GeV/ V10- 39 cm 2 / ' 

(23) 
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M*(M ) 


R*(R Q ) 


Vcsc(V Q ) 




r(yrs.) 


Sun 


1 


1 


4.9 x 10 34 


63 


10 


1.16 


2.49 


8.1 x 10 34 


152 


50 


4.76 


3.24 


6.8 x 10 35 


190 


100 


7.04 


3.77 


1.9 x 10 36 


160 


250 


11.8 


4.60 


6.9 x 10 36 


126 



Table 2: Stellar mass (M*), radius (R*), and surface escape velocity (V esc ) in solar units, 
for metal free stars halfway through hydrogen burning. We have also calculated the capture 
rate C using our fiducial values p x = 10 9 GeV/cm 3 , m x = 100 GeV, and a c = 10~ 39 cm 2 , 
and calculated r using Ca from Table 1. As in Table 1, the entry marked Sun refers to 
the present day Sun and not a zero metalicity star. The capture rate for the Sun still 
uses the fiducial values; in the true present day Sun the true capture rate is much smaller 
(C ~ 10 24 s~ 1 ), mostly due to the much lower DM densities in the solar neighborhood. 



Further, using Eq. (117j) and noting that for the Pop III models of [28j it is roughly true that 
R+ oc M 9,45 , we find that approximately 

C,49x lO^s- 1 (Ml) (-?—) 1 ( E* ^ ( m x \ 1 ( *c \ 

\M Q J Vl0km/s/ V10 9 GeV/cmy VlOOGeV/ Vl0~ 39 cm 2 /' 

(24) 

Table 2 also shows the equilibrium timescale given by Eq.(j3|) using the capture and anni- 
hilation rates determined above. We can see that r is extremely short, compared to the 
life time of a star. Hence for most of the lifetime of the zero metallicity stars, £>t and 
we may use Eq.([5]). 



2.3 Dark Matter Density 

To study the effects of dark matter on the first stars, we need to know the density of the DM 
passing through the stars to determine the capture rate. Simulations have unfortunately 
not (as yet) resolved this issue. Below we will use a variety of DM densities, since these 
numbers are unknown. 



2.3.1 Dark Matter Density before Capture is included 

First we need estimates of the DM density in the region where the star forms, prior to 
including the effects of capture. In a previous paper p] , two of us (with P. Gondolo) used 
adiabatic contraction [15] to obtain estimates of the DM profile. Prior to this contraction, 
we assumed an overdense region of 1O 5 M -1O 6 M with a Navarro-Frenk- White (NFW) 
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profile [16] for both DM and gas, where the gas contribution is 15%. (For comparison, we 
also used a Burkert profile [47] . which has a DM core. A Burkert profile has been shown 
to be a good fit for the dynamics of today's galaxies [4"9| [4"8] ) As the gas collapses, we 
allowed the DM to respond to the changing baryonic gravitational potential, where the gas 
density profiles were taken from simulations of [501 [51] ) . The final DM density profiles were 
computed with adiabatic contraction [M(r)r = constant]. After contraction, we found a 
DM density at the outer edge of the baryonic core of roughly p x ~ 5 (n/cm _3 ) a81 GeV/cm~ 3 
which scales as p x oc r~ L9 outside the core (see Fig. 1 in pQ). Our adiabatically contracted 
NFW profiles match the DM profile obtained numerically in [50J (see their Fig. 2). They 
present their earliest (gas core density n ~ 10 3 cm~ 3 ) and latest (n ~ 10 13 cm~ 3 ) DM 
profiles, as far inward as 5 x 10~ 3 pc and 0.1 pc. The slope of these two curves is the 
same as ours. If one extrapolates them inward to smaller radii, one obtains the same DM 
densities as with our adiabatic contraction approach. The highest DM density found by 
[50] was 10 8 Gev/cm 3 . Should the adiabatic contraction continue all the way to the small 
stellar cores at n ~ 10 22 cm~ 3 (which we doubt), the DM density would be as high as 
10 18 GeV/cm 3 . We note that [52] obtained DM density profiles in galaxies and found that 
adiabatic contraction produces densities that are too high by only a factor of 2 or 3, even 
when radial orbits are included, or in the presence of bars, or in the absence of spherical 
symmetry. Below we will use a variety of DM densities due to the uncertainties. To give a 
sense of the numbers, we compute that in a lOM© dark star, before capture is taken into 
account, there are roughly 10 16 baryons in the star for every WIMP particle. 

2.3.2 Dark Matter Density Including Capture 

The amount of dark matter in the dark star obviously increases significantly due to capture. 
We find that (again for a IOMq star) there are 10 12 baryons for every WIMP particle; i.e., 
the fraction of WIMP particles has grown by a factor of 10 4 due to capture. Of course the 
density of WIMPs is more centrally concentrated, so that the ratio of WIMPs to baryons 
near the center of the star, where the annihilation rate peaks, is higher than these numbers. 
Still, as we will show in the next section, it is remarkable that particles which are 10~ 12 as 
numerous as the baryons can provide the dominant heat source for the star. 

3 Luminosity due to WIMP annihilation: 

We may now compute the luminosity due to WIMP annihilation, 

Ldm = fT A (2m x ), (25) 

where we take the energy per annihilation to be twice the WIMP mass; two WIMPS 
annihilate per annihilation. Here / is the fraction of annihilation energy that goes into the 
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AL(Mq) (Ergs/8) Ldm (Ergs/s) p x (GeV/cm 3 ) for which = Ldm 

Sun 3.9 x 10 33 5.2 x 10 33 7.5 x 10 Y 

10 4.2 x 10 37 3.2 x 10 35 1.3 x 10 10 

50 2.0 x 10 39 2.8 x 10 36 7.5 x 10 10 

100 6.45 x 10 39 7.6 x 10 36 8.3 x 10 10 

250 2.31 x 10 40 2.8 x 10 37 8.5 x 10 10 

Table 3: Luminosity L* due to fusion vs. luminosity Ldm due to DM annihilation in zero 
metallicity stars (as defined in previous tables), using fiducial values for DM properties. 
In the final column we vary p x to determine the value at which Ldm exceeds L*. We 
stress that the DM densities listed here are those in the ambient medium (NFW plus 
adiabatic contraction) rather than the densities after capture; it is these ambient densities 
that determine the capture rate and thus the equilibrium luminosity. 



luminosity. Roughly 1/3 of the annihilation energy is lost to neutrinos that stream right 
out of the star, whereas the other 2/3 goes into electrons, positrons, and photons that are 
trapped in the star and have their energy thermalized. Hence we take [lj 

/ ~ 2/3. (26) 

For t ^> t, which is quickly reached, Eq. (|25p may be rewritten using Eq.(|5]) as 

Ldm = tc(2m x ). (27) 

The WIMP luminosity is given in Table III for a variety of stellar masses together with 
the ordinary fusion-powered stellar luminosity L+ provided by the models of [28J for these 
stars. Roughly, using Eq. (|24p we may write 

L DM = 5.2 x 10 33 erg/s ( ^ ( ^ *) ( %-^) ( ' . (28) 

1 \M Q ) Vl0 9 GeV/cm 3 y VlO^cm 2 / V ; 

The WIMP luminosity depends linearly on the the WIMP density passing through the stars. 
We have also computed the WIMP energy density p x , C rit that is required in order for the 
WIMP annihilation energy to equal the ordinary stellar luminosity, which will dramatically 
alter the properties of the first stars [29J. We stress that the DM densities listed here are 
those in the ambient medium (NFW plus adiabatic contraction) rather than the densities 
after capture; it is this ambient density that determines the capture rate, and also therefore 
the the annihilation rate (in equilibrium) , and consequently the luminosity. For any WIMP 
densities higher than this value, the star's luminosity is dominated by annihilation energy 
(rather than by ordinary fusion). We note that the luminosity due to fusion in the zero- 
metallicity stars [28] scales as L* oc M 2 (see figure 1) whereas Ldm °t Af , so that DM 
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heating is relatively more important in lower-mass stars. The DM luminosity is dominant 
for the values of /o x , C rit given in the last column of Table III, and is dominant for all relevant 
stellar masses for /) x >9x 10 11 GeV/cm 3 ; thus for dark matter densities in excess of this 
value, it will be the DM heating that determines stellar properties, and the star is a dark 
star. The two luminosities are plotted in Figure 1 where they can be compared with one 
another. If the first stars are observed (e.g., by James Webb Space Telescope) to have the 
properties predicted of fusion-driven stars, then one could use these results to constrain 
WIMP properties to make sure that DM annihilation remains subdominant. 

A potentially important consideration is the time for which a dark star will be able to burn 
DM at the calculated rate, given that it will annihilate its captured store of DM in a time 
of order r as given by Eq. . The maximal burning time is determined by the total mass 
of DM in the region of phase space that intersects the star. This mass could be small if 
the star were fixed exactly at the central cusp of the DM halo. However, we expect that 
given the complexities of the collapse process, the star will have some nonzero velocity and 
thus 'wander' through a region of some radius significantly larger than that of the star. 
This makes much more DM mass available for burning, at the expense of somewhat lower 
average DM density. (It should then be noted that the fiducial DM density numbers we 
assume are not those expected in the innermost core but over some significantly larger 
region.) Estimation of the expected degree of wandering, the detailed DM density profile, 
and therefore the maximal burning time of the dark stars will be left for future and more 
detailed study. 



4 Eddington luminosity 



We as yet know little about the effect that significant dark matter heating would have on 
the structure of zero-metallicity stars, but we can nonetheless place an approximate upper 
bound on the stars' mass if we assume that they must be sub-Eddington. In the current 
paper, we simply take the ZAMS stars, add the DM luminosity, and ask whether or not 
the resulting dark stars are self-consistent. In other words, is it possible for dark stars 
of this mass to exist? We compute the Eddington luminosity L^d for these objects, and 
ascertain whether or not the DM luminosity is in excess of this value. If Ldm > -^Edd; then 
the star of this mass cannot exist: the pressure from the star would be so large as to blow 
off some of the mass. In reality one should do a different problem: one should really follow 
the protostars and compute their structure as they accrete mass on their way to becoming 
ZAMSs. Such a calculation would essentially compare accretion luminosity (the value of 
which would depend on the nature of the accretion, e.g. spherical or disk accretion) to 
the Eddington luminosity; here we are comparing dark star luminosity to the Eddington 
luminosity. 
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The Eddington luminosity is defined (e.g. [H]) to be 



L E dd = *, (29) 

where G is Newton's Constant, c is the speed of light, M* is the mass of the star, and k p 
is the opacity of stellar atmosphere. Since the first stars' stellar atmospheres are hot and 
nearly metal-free, the opacity is dominated by Thompson scattering, so we take 

L Ed d = 1-4 x 10 38 erg/s {MJM Q ) = 3.5 x 1O 4 (M*/M )L . (30) 



In Figure 1, we have plotted three luminosities as a function of stellar mass for zero metal- 
licity stars: the luminosity due to fusion, the luminosity due to WIMP annihilation (for a 
variety of WIMP densities in the star) , and the Eddington luminosity. Since the Eddington 
luminosity scales as L oc M* whereas the DM luminosity scales as Ldm ck M*' 55 , for a large 
enough dark matter density the two curves will cross for some stellar mass. The lightest 
stellar mass for which Ldm > -^Edd then constitutes an approximate upper mass limit to 
the first stars, because the star of that mass is unable to accrete any further due to the 
radiation pressure from the WIMP annihilation. If one assumes that accretion efficiently 
drives up the mass of any Pop III star, then the mass of the first star will be uniquely 
determined by the properties of the DM. Using Eqs. (f28l) and (f29|) . we get: 



For example, we find that for a dark matter density of p x = 1.1 x 10 13 GeV/cm 3 , the first 
stars cannot be more massive than 1M . Figure 2 illustrates M™ ax as a function of WIMP 
energy density for several values of the scattering cross section. 

In principle, these arguments can be turned around to place a bound on the scattering cross 
section. Figure 3 is a plot of WIMP scattering cross section vs. WIMP mass. Experimental 
bounds are shown. The horizontal lines indicate the values of a c that correspond to different 
values of p x at which a 1 M star is Eddington limited by DM annihilation. (The scaling 
to other stellar masses is straightforward since Ldm cx .) Once the mass of the first 
stars is determined, then one could rule out any combination of a c and p x for a given WIMP 
mass that would preclude stars of such a mass from forming. As an extreme example, if 
one were to believe adiabatic contraction all the way to the limit where the protostellar 
core has gas density n ~ 10 22 cm -3 , the WIMP density would reach p x = 10 18 GeV/cm 3 . 
In this case the Eddington limit would be reached for masses <C 1 M for our fiducial 
scattering cross section. If, instead, the first stars are observed to form with masses larger 
than one solar mass, and such an enormous WIMP density were found to be sensible, one 
could place a bound of a c < 3 x 10 4 cm 2 on WIMPs for almost any mass (the equilibrium 
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DM luminosity is independent of m x ), which would be the tightest known bound on the 
spin dependent scattering cross section by several orders of magnitude; see figure 3. In this 
extreme case we could also put interesting limits on spin independent scattering; see figure 
4. 

Obtaining bounds on the WIMP parameters in this way depends on detailed understanding 
of the ambient WIMP density within which the dark star resides. Clearly this is a very 
difficult problem which will not be solved in the near future. Indeed, it is likely that 
the 1O 6 M0 haloes containing the dark stars will merge with other haloes and that the 
stars will not remain forever in regions of high DM density. However, simulations of 
structure formation are becoming ever better and one may hope someday to address this 
question. 

5 Conclusion 

In summary, Pop. Ill stars are expected to reside at the core of a 10 5 - 10 6 M Q dark matter 
halo. Previously two of us pQ discussed the importance of dark matter annihilation in the 
first stars, proposing the existence of "dark stars" powered by DM heating. Even once 
the DM initially collapsing with the baryons into the star runs out as a source of fuel, it 
can be replenished by capture of more DM from the ambient medium (as the DM passes 
through the dark star). In this paper, we have estimated the rates of WIMP capture and 
self-annihilation in such stars. We have shown that when these rates are in equilibrium, 
the accompanying heating would provide an energy source that can rival nuclear fusion 
and prolong the dark star phase, if the core dark matter density is sufficiently high. Such 
densities seem plausible based on analytic models of adiabatically contracted halos [I], 
suggesting that DM heating may radically affect the structure of the first stars. 

The two key uncertainties in this work are: (i) the scattering cross section must be at (or 
near) the experimental limit, and (ii) the ambient DM density must be high enough for 
capture to take place. The second criterion is likely to be true for a while after the star is 
created, but it is not clear for how long. Once the halo containing the Dark Star merges 
with other objects, it is not clear how long the central DM in the halo remains undisturbed, 
and it is not clear how long the Dark Star remains at this central point. In principle the 
capture could continue indefinitely so that a dark star could still exist today, but this is 
very unclear. 

For high enough DM density, DM heating will lower the Eddington stellar mass limit to 
provide an upper mass cutoff for Pop. Ill; because DM heating might also affect the 
formation properties of Pop. Ill stars, the gross properties of these objects may well be 
determined by the particle properties of dark matter. And conversely, inferred properties 



15 



of Pop. Ill stars (or even future direct observations) might be used to strongly constrain 
DM masses and interaction cross sections. 
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Log[M,/M J 



Figure 1: Log-Log plot of stellar luminosity vs. stellar mass (in solar units). The dashed 
line shows Eddington luminosity for Thompson scattering. Triangles show the data points 
for luminosity of zero metallically stars in Table 3; the solid line is a fit to these points. The 
remaining lines indicate the DM annihilation luminosity for different DM density assump- 
tions. We stress that the DM densities listed here are those in the ambient medium (NFW 
plus adiabatic contraction) rather than the densities after capture; it is these densities that 
determine the capture rate and the luminosity. The dotted line is our fiducial example with 
p x = 10 9 GeV/cm 3 and a c = 10 -39 cm 2 (squares are from the third column of Table 3.) If 
we keep the cross section fixed, the middle dot-dashed line corresponds to a DM density 
of 10 12 GeV/cm 3 , and the top diamond line would have a DM density of 10 14 GeV/cm 3 ; 
in this final case, the DM luminosity would dominate over the Eddington luminosity for a 
star with a mass of < 1 M & . 
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Figure 2: Maximum stellar mass due to the Eddington limit as a function of p x for a fixed 
cross section. The different lines correspond to different spin-dependent cross sections: 
a c = lCT 39 cm 2 (solid line), a c = 10~ 41 cm 2 (dotted line), and a c = 1CT 43 cm 2 (dashed 
line). 
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Figure 3: Bounds on spin-dependent scattering cross section as a function of WIMP mass 
from experiments, as well as potential bounds from the effect of WIMP annihilation on 
the first stars. Displayed limits include various direct and indirect detection experiments 
as labeled; the tightest present bounds on spin-dependent cross-sections are from Super- 
Kamiokande (labeled SuperK), at around 10 -38 cm 2 [311 [32], [351 ESI EI] . The three horizonal 
lines are bounds that would result from the discovery of Pop III stars of 1M@ for a variety 
of DM densities inside these stars (due to the fact that DM annihilation pressure would 
otherwise prevent their formation). The lines are labeled as follows: DS I corresponds to 
ambient density p x = 1.1 x 10 13 GeV/cm 3 and would lead to the bound a c < 10~ 39 cm 2 ; 
DS II corresponds to p x = 3 x 10 15 GeV/cm 3 and would lead to a c < 10 -41 cm 2 ; and DS 
III corresponds to p x = 10 18 GeV/cm 3 and would lead to a c < 3 x 10~ 44 cm 2 . 
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Figure 4: Same as previous figure but for WIMPs with spin-independent interactions 
[381 EH [341 [37]. 




22 



